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Abstract. We give explicit, polynomial-time computable formulas for the number of integer points in any two- 
dimensional rational polygon. A rational polygon is one whose vertices have rational coordinates. We find that the 
basic building blocks of our formulas are Dedekind-Rademacher sums, which are polynomial-time computable finite 
Fourier series. As a by-product we rederive a reciprocity law for these sums due to Gessel, which generalizes the 
reciprocity law for the classical Dedekind sums. In addition, our approach shows that Gessel's reciprocity law is a 
special case of the one for Dedekind-Rademacher sums, due to Rademacher. 



The full beauty of the subject of generating functions emerges only from tuning in on both chan- 
nels: the discrete and the continuous. Herb Wilf |WI p. vii] 



1 Introduction 

We define a two-dimensional polytope V as a compact subset of M 2 bounded by a simple, closed 
polygonal curve. V is called a rational polytope if all of its vertices have rational coordinates. We 
give explicit, polynomial-time computable (in the logarithm of the coordinates of the vertices) 
formulas for the number of integer points in any two-dimensional rational polytope and its integral 
dilations. We emphasize an expository flavor in this paper. In the current literature, there are 
either 'nice' formulas that do not appear to be polynomial-time computable |BVl lDR KKl IKPl IP] . 
or there are polynomial-time computable algorithms without 'nice' formulas |Baj . Asking for both 
seems to be asking for too much; but in M 2 we show that we can have our cake and eat it, too. 

To fix notation, let V° be the interior of V, and V = V be the closure of V . For t € N, let 
L(P°,t) = # (tV° fl Z 2 ) and L(V, t) = # (tV n I?) be the number of lattice points in the interior 
and closure, respectively, of the dilated polytope tV = {{tx, ty) : (x, y) E V}. Ehrhart, who initiated 
the study of the lattice point count in dilated polytopes jEj, proved that L(V°,t) and L(V,t) are 
quasipolynomials in t. A quasipolynomial is an expression of the form c n (t) t n + • • • + c\(t) t + co(i), 
where cq, . . . , c n are periodic functions in t. 

A natural first step is to fix a triangulation of V, which reduces our problem to counting integer 
points in rational triangles. However, this procedure merits some remarks. First, triangulation 
is not easy, but there has been remarkable progress recently, so that we can triangulate a convex 
polygon with n vertices in roughly n steps On the other hand, we are concerned with the 
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efficiency of our formuias with respect to the coordinates of the vertices of V, and not the number 
of vertices. 

Another point of non-trivial significance is the number of lattice points on rational line segments, 
namely the boundary of our triangles. Although this is considerable easier than enumerating lattice 
points in 2-dimensional regions, it is still non-trivial and has only recently been completely solved 
BR, T . It is amusing that counting lattice points on line segments gave rise to links with the 
Frobenius coin-exchange problem and the number of representations of an integer by a linear form. 

After triangulating V, we can further simplify the picture by embedding an arbitrary rational 
triangle in a rational rectangle: 



Since rectangles are easy to deal with, the problem reduces to finding a formula for a right- 
angled rational triangle. Such a formula is given in Section [2] using generating functions; this 
derivation is a refinement of a previously introduced method |Bej . We find (Section |2J that the 
basic building blocks of the lattice point count formulas for any two-dimensional rational polytope 
are the sawtooth function ^ 

((»)) :=x-[x]-- 

and the Dedekind-Rademacher sum 



k=0 



Here a and b are integers, and t is a real number. We use the name Dedekind-Rademacher sum 
in a somewhat lenient fashion; often {{x)) is defined to be if a; € Z, also Rademacher's original 
definition is |Di| IMe| IRaj 



( ( a(k + y) ^\\ ( ( k + y 



s(a,b;x,y) := ^ ( ( hx 

k=0 

Here a and b are integers, whereas x and y are real. However, it is clear that the different use 
of ((. . . )) only results in a difference of the arithmetic sums in a trivial term. Also, a and s are 
strongly linked via 

a(a, b,t) = s ( a, b; — , 
V « 

and 

y 

s(a, b; x, y) = a (a, b, ay + bx) + ^((ya + xb)) . 

We chose to use a rather than s because of its natural appearance in our formulas. There exists a 
two-term reciprocity law for these sums (KllRaj . which enables us to compute a(a, b, t) in polynomial 
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time, similar in spirit to the Euclidean algorithm. From this fact we conclude that our lattice point 
enumerator for V is polynomial-time computable (Section QJ. 

As a by-product of our formulas, we rederive in Section |S] two reciprocity laws as corollaries: the 
two-term law for the classical Dedekind sum ( |Dej . Chapter 2 of |RG| ). and a two-term law for 
generalized Dedekind sums due to Gessel |U| . In fact, our approach shows that Gessel's reciprocity 
law is a special case of the reciprocity law for Dedekind-Rademacher sums, a theorem due to 
Rademacher |R,aj . 



2 Generating functions 



In |Be| . the first author used the residue theorem to count lattice points in certain tetrahedra. Here 
we adjust and expand these methods to the rectangular triangles we reduced the discussion to in 
the introduction. Such a rectangular triangle T is given as a subset of M. 2 consisting of all points 
(x, y) satisfying 

a b 
x > ^, y > -, ex + jy <r 

for some integers a, b, d, e, /, r with ea + fb < rd. Because the lattice point count is invariant under 
horizontal and vertical integer translation and under flipping about x- or y-axis, we may assume 
that o, b, d, e, /, r > and a,b < d. Let's further factor out the greatest common divisor c of e and 
/, so that e = cp and / = cq, where p and q are relatively prime. Hence 



(2) 



T= |(x,y) G R 2 : x > — , y > —, cpx + cqy < r 
To derive a formula for L (T,t) we interpret, similarly as in |Bej . 



L(T,t)=#Um,n)e 



ta 



tb 



■ Tn ^ — - , n > — , cpm + cqn < tr 
d d 



as the Taylor coefficient of z tr of the function 

( \ I 



E 



cpm 



m> 



E 



.,cqn 



+1 




1 - Z C P 



I- z cq l_ z 



y U + V 



(1 - z c p) (1 - z c i) (1 - z) ' 
where we introduced, for ease of notation, 

'ta - 1" 



u :- 



d 



+ l)cp 



and 



v := 



tb- 1 

d 



+ l)cq. 



(3) 



(4) 



We present two methods on how to extract the lattice point count from this generating function: 
partial fractions and the residue theorem. Both are inspired by works on generalized Dedekind 
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sums, the first one by Gessel [U], the latter one by Zagier [Zj. In fact, both ways are completely 
equivalent, since our generating function is rational. However, to please both algebraically and 
analytically minded readers, we give two proofs of the following 

Proposition 2.1 For the rectangular rational triangle T given by 

L{T,t) =-^(tr-u-v) 2 + l(tr-u-v) (± + ± + J-\ 
2&pq 2 \cp cq c A pq J 

1/ 1 1 \ 1 (p q 1 

+ T 1 + — + — - + - + 



4 \ cp cq J 12 \q p c 2 pq 

1 1 u + v-tr\ \- tr 1 v-^ \- tr+1 



1 v ^ A^" tr 1 x ^ X u ~ tr 

+ ~c~p~ ^ (1 - \ c i) (1 - A) + ~q ^ (1-A c p)(1-A) ' 

where u and v are given by 

It will be useful to have the Laurent expansion of the factors of our generating function. The 
following lemma will provide a bridge between the residue method and the partial fraction method. 



Lemma 2.2 Let a,b be positive integers, and A° = 1. Then 

( z ~ A ) + -TT-r- + 0{z - A) . 



1 - z ab ab v ' v ' 2ab 
Proof. First, 

Res( r, z = A ) = lim r = • 

V 1 - z ab ) 1 - z ab ab 

For ab = 1, the statement is trivial, so we may assume ab > 2. Then the constant term of the 
Laurent series of ^ a b can De computed as 

1 A \ ab(z - A) + A (l - z ab ) ab - ab\z ab - 1 

hm -t + — — = hm — — -t--^- = hm 



A VI - z ab ab(z - A) J ab(z - A) (1 - z ab ) z->\ ab (1 - z ab - (z - X)abz ab - 1 ) 

-X(ab - l)z ab - 2 ab-1 
z->\ -2abz ab - x -(z- \)ab{ab - l)z ab - 2 2ab 

□ 

First proof of Proposition (partial fractions). To make life easier, we translate the coefficient 
of z tr of our generating function, which yields the lattice point count, to the constant coefficient of 
the function 



(1 - z°p) (1 - z c i) (1 - z) 



(5) 
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This is a proper rational function because 



a b 
cp — h cq- < r 
a a 



(T ^ 0!), which implies 
u + v — tr — cp — cq — 1 



ta - 1 
d 



cp 



tb-1 
d 



ta tb 

cq — tr — 1 < — cp H — -cq — tr — 1 < — 1 . 
d d 



By expanding (JSJ) into partial fractions 

„u+v—tr 



(1 - 2°P) (1 - 2^) (1 - 2) 



E 



+ E A+ E 



z — A ^— ' z — A 



£>A 



_C\ 

z - a ^ (z - xy 



3 j-, tr—u—v _ 



fe=l 



fe=l 



we can compute L (T, i) as the constant coefficient of the right-hand side: 



L(T,t) 



E t - E t + E + 3 



Ca -Da 
A A 2 



(6) 



The computation of the coefficients A\ for X cp = 1 ^ A c is straightforward: 



yl\ = lim 



(z - X)z u+v - tr 



A 



v—tr 



lim 



(z-A) 



A 



u— tr+1 



♦A (1 - z cp ) (1 - z c i) (1 - z) (1 - A c <?) (1 - A) z— >A 1 - z c p cp (1 - A c <?) (1 - A) 
Similarly, we obtain for the cgth roots of unity X cq = 1 7^ A c 

^m— tr+l 



B> 



eg (1 - A c p) (1 - A) 



The coefficients D\ and C\ are the two leading coefficients of the Laurent series of © about a 
nontrivial cth root of unity A. Using Lemma 12.21 they are easily seen to be 



X 



-tr+2 



c 2 pq(l — A) 



and 



° X ~ { ~2c~p " 2c~q + 



u + v — tr + 1\ X 



c 2 pq 



i-tr+1 



+ 



A 



-tr+2 



1 - A c 2 pq(l - A) 2 ' 

Finally, we obtain the coefficients E\~ from the Laurent series of © about z = 1 (by hand or, 
preferably, using a computer algebra system) as 



E-i 



c 2 pq 



u + v — tr + 1 1 1 

e 2 = 2 h o — h 

c^pq Zcp Icq 



and 



Eh 



(u + v — tr) 2 u + v — tr 



2c 2 pq 



+ 



1 



1 1 



+ — + 



+ 



1 ( 1 



c 2 pq cp cq) 4 \cp cq 



+ 



1 



1 I p 1 q 



12 \q c 2 pq p / 
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Putting these ingredients into © yields the statement. □ 

Second proof of Proposition \2.1\ (residue theorem). The sought-after Taylor coefficient of © can 
be shifted to a residue: 

( u+v—tr—l \ 

(i-*»)(i-z-)(i-z) -' = °) (7) 

If the right-hand side of (JJJ) counts the lattice points in tT, then what we have to do is compute 
the other residues of 



/(*) 



gU+v— tr— 1 



(1 - Z C 'P) (1 - z c i) (l-z) ' 

and use the residue theorem for the sphere C U {oo}. Aside from 0, / has poles at all cpth and cgth 
roots of unity; note that the nonemptyness of T implies Res(/(z), z = oo) = 0. 

The residue at z = 1 can be easily calculated as 

Res(7(»,z = l) = Res(e z f(e z ),z = o) = + « - trf 

1 , . / 1 1 1 \ 1 / 1 1 \ 1 fp q 1 

-f- (« + « - tr) — + — + -=— - 7 1 + — + — -T7T - + - + 



2 \cj) q c 2 pqj 4 \ cp eg/ 12 \q p c 2 pq 

It remains to compute the residues at the nontrivial roots of unity. Let A be a nontrivial cth roots 
of unity. Putting the Laurent expansions of the different factors of / together, the residue of / at 
A can be easily derived via Lemma 12.21 as 

„ /, , , x \ / u + v-tr 1 1 \ \~ tr \- tr+1 
Res(f(z),z = X) = = + _ + _ _ _ + 



c 2 pq 2cp 2cq J 1 — A c 2 pq(l — A) 2 

Note that we have to add this expression over all nontrivial cth roots of unity. 

Now let X cp = 1 ^ X c . Since p and q are relatively prime, / has a simple pole at A, whose residue 
can be determined easily using Lemma 12.21 

Res( f{z),z = A) = — — — Res —, Z = X 



(1-A C 9)(1-A) yi-Xcp' J cp(l - A C( ?) (1 - A) ' 

Again, we have to add up all these A's. Finally, we obtain a similar expression for the cgth roots 
of unity, and the statement of the proposition follows by rewriting by means of the residue 
theorem. □ 



In the following section, we will further describe the finite sums appearing in the lattice point count 
operators; consequently, we will be able to make statements about their computational complexity. 
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3 Using the Dedekind— Rademacher sums as building blocks 

We will now take a closer look at the finite sums over roots of unity appearing in Proposition 12.11 
namely, 

1 x \ w 



cp V, (1 - xcq ) (1 - A ) 

for some integers c,p, q, w, where p and q are relatively prime. The fact that this is a finite Fourier 
series in w and the appearance of two factors in the denominator suggest the use of the well-known 
Convolution Theorem for finite Fourier series: 



Theorem 3.1 Let f(t) = — ^ a A A* and g(t) = — ^ b\X t . Then 

x N =i \ N =i 



A Ar =l m=0 

□ 



The key ingredient to be able to apply this theorem to our case is 
Lemma 3.2 For p,t G Z, 



L(J,t) 



+ 1. (8) 



p ^— ' A — 1 \ \ P J J 2p 

Recall that {{x)) = x — [x] — 1/2. This lemma is well-known (see, e.g., ,RG , p. 14), however, for 
sake of completeness we give a short proof based on the methods of section 

Proof. Consider the interval 2 := [0, -], viewed as a one-dimensional polytope. Then the lattice 
point count in the dilated interval is clearly 

~t~ 

_p_ 

On the other hand, we can write this number, by applying the ideas in section as 

L(Xt)= Res ( (1 _" P)(1 _ z) ^ = / 

(Equivalently, we could expand this generating function into partial fractions.) Using again the 
residue theorem, this can be rewritten as 

/- \ t 1 1 1 A"' 
L ^=p + Yp + 2-p £ A^T- < 9 > 

Comparing (JSJ) with @ yields the statement. □ 
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Corollary 3.3 For c,p, q,t GZ, (p, q) = 1, 

1 E - 

Here, qq -1 = 1 mod p. 



i v _* J-((=S0)-± *<* 

C P \c P= i-£\c * q [ otherwise. 



Proof. If c|t, write i = cw to obtain 

1 A* 1 ^ X cw 1 ^ X w 1 A 9 " 1 " 1 

cp ^ 1 - X c i ~ cp ^ 1 - \ c i ~ p ^ 1 - \i ~ p ^ 1 - A 

^ A c p=1^A c ^ A c p=1^A c AP=1^A ^ AP=1^A 

(*) / /" — g -1 ^ \ \ 1 / / -g _1 A\ 1 



p J J 2p \ \ cp J J 2p 
Here, (*) follows from Lemma 13.21 If c does not divide t, let £ = e 2m ^ cp . Then 

P- 1 c- 1 ^{mc+np)t 1 C_1 , t £ mc * 



i \t i £i ^ ^ c( mc + n P)t i r *: 

_L y A = _L y i_ = — V £ npt V - 

co Z- / ^ _ \cq CD Z_^ Z-/ 1 _ f(mc+np)cq cv Z-^ ^ Z-/ 1 _ tmc 2 q 
y A C P=1^A C ' r m=lti=0 s ^ n=0 m=l s 



. 



□ 



Corollary 3.4 For c,p, g,i £ Z, (p, q) = 1, 

* A J^. (! " AC<? ) (1 " A ) = ~° c) " ((^)) + ^ ((c 

Proof. We will repeatedly use the periodicity of the sawtooth function. One consequence is, for 

g((=^))-«-». 

the proof of which is left as an exercise ( [KG] , p. 4). Now by Lemma 13.21 
1 V "A A^ 1 V A^ 1 V A^ 

cp ^ (1 - A) ~ cp ^ (1 - A) ~ cp ^ (1 - A) 

^)) " 2c~p ~ p (" ((t)) " 2c) = ~ ((^)) + p ((t)) ' 
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Finally we use the Convolution Theorem 13. II and Corollary 13.31 to obtain 



- y - 



<* 0- ~ A-) (1 - A) 

Cp— 1 / , , — 1\\ -i \ / / / I \ \ \ i / / / \ 

' II- x m \ \ I \ / / / —(t — m } \ \ I I I —[I — m j 



m = 

elm 



S( (( C p )) 2p) ( (v C P J/P 



- — " / / -q 1 k\\ / / k t \\ 1 1 / / — g \ f f —t 



P / J \ \P C P J J V — \\ P / / V V c 
fc=0 fc=o 

p_1 n\ i t^ff-t 



+ O n ^[[ r, ml) 7r,2^ 



2pf^\\p cp J J 2p^ o \\c 
= £ ((~P~)) ((7 " cp")) + 2~P ((")) + 2p~ ((")) ~ 2p~ ((~ 

~^((p)) ((~P "<?>)) " (("cp)) + 2p~ ((t 

The statement follows now by definition of the Dedekind-Rademacher sum ([T|). □ 

This corollary describes the finite sums in Proposition 12.11 One of them actually turns out to be 
of an even simpler form. To show this, we first need to rewrite Proposition 12. II for the special case 
where T has the origin as a vertex: 

Proposition 3.5 For the rectangular rational triangle T given by PJ) with a = b = 0, c = r = 1, 

and p and q relatively prime, 

T _ , t 2 t ( 1 1 \\ 1 1 (p q 1 
L (T,t) = — + - - + - + — +- + — - + - + — 
ipq 2 \p q pqj 4 12 \q p pq 



-a(q,p,t) -a(p,q,t) 



3 - G 



Proof. Proposition 12.11 gives for this special case 

L(T,t) = — + -(- + - + -) + + - + -)+ -(- + - + - 

Ipq 2 \p q pq ) 4 \ p g/ 12 \g p pg 

+ P ^ (1-A«)(1-A) + g ^ MM • 
The statement follows now with Corollary 13.41 □ 
We use this Proposition to show 
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Lemma 3.6 For p,t G Z, 



p-i 

k=0 



k + t 



.p_+L+p((i 

24 6p 2\\p 



Proof. Consider the triangle A := {(x,y) € 



+ py < l} and its integer dilates. By summing 



over vertical line segments in the triangle, we obtain 



L (A, i) = (t - pm + 1) = (t + 1) 

m=0 

t 2 /l l\ 1 p 

2p \p 2J 2 8 



+ 1 



P 
2 



KG 



On the other hand, we can compute the same number via Proposition 13.51 



x * * / 2 , \ 1 1 

L(A,i= — + - - + 1+- + — lp + 
2p 2 \p ) 4 12 \ p 



cr(l,p,t) 



1 



+ 1 



+ 



Again we used pi)|). Equating (jTTj) with ((T2*]) yields the statement. 

Using these ingredients, we can finally restate Proposition 12 . 1 1 as our main theorem: 



(11) 



(12) 



□ 



Theorem 3.7 For the rectangular rational triangle T given by |]|), 



L{T,t) 



-g— (tr - u - vf + (tr - « - v) — + — + -g— + - 

2c z pq \ 2cp 2cq c z pq cpq \ \ c J 



1 



1 



1 



1 



+! + ![£ + £, , , , . 

4 12 \q p) 2Apq c 2 pq \2p 2q 



1 + > + fl4)^ 



tr — v 



cp 



tr — u 
cq 



+ - 



1 



tr - 1 



.Here it and v are given by J^p. 



+ 



2pg 



+ 
tr - 1 



1 1 1 

cpq 2p 2q 

2 



tr — v 



o- q,P, 



a p 



tr\ 



tr — u 



Proof. By Lemma 13.21 



1 x ^ X w 
c ^ 1-A 



By Corollary 13.41 and Lemma 13.61 



- y - 

c . (1 



c " (1-A) 2 



-W 
C 



-a (1, c — w) 



1 

2^ 



c 



1 

4^ 



(13) 
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Now simplify the identity in Proposition 12. II bv means of (|13|). (|14|) . and Corollary 13.41 □ 

We mention the work of Brion and Vergne [BV on rational polytopes, where they give formulas 
for the lattice point enumerator of a convex rational polytope in terms of certain Todd differential 
operators. These are interesting connections to topology. The salient difference in approaches lies 
in our explicit use of Dedekind-Rademacher sums and their computational efficiency. Although 
|BVj offer an approach which might eventually yield similar results, the analysis of the equivalence 
between our different approaches would be the content of another paper. The results in |BV| main 
result on p. 831] are coordinate-free, and have their theoretical merits. In contrast our results are 
coordinatized and are useful in the context of computer-related applications. 

4 Remarks and consequences 

Rademacher's original definition |Raj of his generalization of the Dedekind sum is 

— 4g((^H)'(m)"' 

defined for a, b € Z, x, y € R. Here, 

(M) *.-! ((*)) if*^ 
[{X>> — \ if x <E Z 

is the sawtooth function which also appears in the classical Dedekind sum S(a, b\ 0, 0). The impact 
of the slightly different definition of the sawtooth function is not crucial for our lattice point count 
formula. In fact, it is easy to see that 

<r(a,M)=s(a,fc;~,o) " 5 ((£)) ■ (IB) 

An important property of S(a, b; x, y) is Rademacher's reciprocity law |Ra| 

S(a, b; x, y) + 5(6, a; y, x) = 

= J -j + A(f + 3 + 1) ifbothx,y G Z 
\ (Oz))*((y))* + \ HMV) + ITbMay + bx) + %ih(xj) otherwise, 

where ^ 

ifj 2 (x) := (x - [x]) 2 - (a? - [a?]) + - 

D 

is the periodic second Bernoulli polynomial. The equivalent reciprocity law for a (a, b, t) was first 
presented in [Kj we will rediscover it in the last section. 

Among other things, Rademacher's reciprocity law allows us to compute S(a, b; x, y) (and hence 
ct(o, b, t), the nontrivial part of our lattice point count formulas) in polynomial time, by means of a 
Euclidean algorithm using the first two variables: simply note that we can replace a in S(a, b; x, y) 
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by the least residue of a modulo b. It is amusing to note that a(a, b, t) appears in the multiplier 
system of a weight-0 modular form |Roj . 

To complete the picture for an arbitrary two-dimensional rational polytope V, we return to the 
statements in the introduction. After triangulating V, the problem reduces to rational rectangles 
and the rectangular triangles which were treated above. A lattice point count formula for a rational 
rectangle K is easy to obtain: suppose K has vertices (^, Sf) , , , (^f, %f\, with 

aj < bj, then it is not hard to see that 



( 


~tbi 




ta\ — 1 


)( 


~tb 2 ~ 




ta 2 - 1 




d _ 




d 




d _ 




d 



Together with the above established remark on computability of a(a,b,t), we summarize some 
statements in 



Theorem 4.1 Let V be a two-dimensional rational polytope. The coefficients of L(V,t) can be 
written in terms of the sawtooth function ((...)) and the Dedekind-Rademacher sum a(a,b,t). 
Consequently, the formula given by Theorem \'J. l\ for the lattice point count operator can be computed 
in polynomial time in the logarithm of the denominators of the vertices ofV. □ 

Barvinok |Baj showed that for any fixed dimension the lattice point enumerator of a rational 
polytope can be computed in polynomial time. The distinction here is that we get a simple formula, 
which happens to be polynomial-time computable. 

In dimensions greater than 2, things get more complicated. We can still get formulas through the 
methods introduced here; however, even the existence of (possible) three-term reciprocity laws for 
functions appearing in the lattice point count does not guarantee polynomial-time computability. 
The details will be described in a forthcoming paper BDR . 



5 Reciprocity laws 



As another remark, we can recover the reciprocity law for the classical Dedekind sum ( |Dej . Chapter 
2 of RG ) from our formulas: 



Corollary 5.1 (Dedekind) 



s( , 6; o,o) +S (M;0,o) = -i + i(^ + i + ^ 



Proof. It is well known EJ that the constant term in the Ehrhart polynomial, the integer-point 
enumeration function of a lattice polytope (that is, a polytope with integral vertices) equals the Euler 
characteristic of the polytope. Consider the simplest case of our triangle mentioned in Proposition 
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13.51 If we dilate this polytope by t = pqw, that is, only by multiples of pq, we obtain the dilates of 
a lattice polytope V. Proposition 13.51 simplifies for these t to 

L(V,w) =P^ + ^(p + q + l) + ] + !-(?- + + -a(q,p,0)-a(p,q,0)+l . 
v / 2 2 4 12 \q p pq J ' ' 

On the other hand, we know that the constant term is the Euler characteristic of V and hence 
equals 1, which yields the identity 

7 + 7^ (- + - + — ) ~cr(q,p,0) -a(p,q,0) = . 
4 12 \q p pq J v 

The statement follows now by rewriting the Dedekind sums in terms of the original definition via 

csd. □ 

As a concluding consequence of our formulas, we rederive a reciprocity law due to Gessel 
interpreting it at the same time geometrically. 

Corollary 5.2 (Gessel) Let p and q be relatively prime and suppose that t is an integer such that 
1 < t < p + q. Then 

1 A* 1 A* 

+ 



t 2 t (\ 1 \\ 1/1 1 \ 1 (p I 



2pg + 2 Vp + g + pgj 4\p + q +1 ) 12 \q + pq + p 

It is easy to see that the reciprocity law for classical Dedekind sums (Corollary 15.11) is a special case 
of Gessel's theorem. Before proving Gessel's theorem below, we find it useful to have the lattice 
point count operator for the interior of our polytope. The following central theorem, conjectured 
by Ehrhart [E] and first proved by Macdonald |Maj . enables us to do this: 



Theorem 5.3 (Ehrhart— Macdonald Reciprocity Law) If a rational polytope V is homeomor- 
phic to a d-sphere then 

L{V°-t) = {-l) d L{V,t) . 

□ 

Note that this theorem allows us to conclude a computability statement for the lattice point count 
in the interior of a two-dimensional rational polytope similar to Theorem 14.11 Using Theorem 15.31 
we get from Proposition 12. II the 

Corollary 5.4 For the rectangular rational triangle T given by 0j with a = b = 0,c = r = l, and 

p and q relatively prime, 

_ frr0 , t 2 t (I 1 1 \ 1 / 1 1\ 1 (p q 1 



+ 







1 






V 












\P 




pqj 





2pq 2 \p q pq J 4\ p q) 12 \q p pq 
A* 1 ^ A* 



l ^ A 1 

v ^ (1-A9H1-A) + a ^ IT 
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□ 

Alternatively, we could have derived Corollary 15 .41 from scratch in a similar way as Proposition 

Proof of Corollary \5.2L Consider dilates of the triangle given in Corollary 15. 4[ that is, 

tT° = { (x, y) € M 2 : x, y > 0, px + qy < t} . 

By the very definition, tT° does not contain any integer points for 1 < t < p + q, in other words, 
L (T°,t) = 0. Hence Corollary 15.41 yields an identity for these values of t: 

t 2 t (\ 1 1 \ 1 / 1 1\ 1 (p q 1 
= - =(- + - + — +-r l + - + - +T^ - + - + — 



2pq 2 \p q pq J 4 \ p q) 12 \q p pq 
1 x ^ A* \ x ^ \ l 

+ p 2^, (l-\i)(l-X) + q^, (1 - AP) (1 - A) ' 

We can rephrase Corollary 15.21 in terms of Dedekind-Rademacher sums by means of Corollary 13.41 



Corollary 5.5 Let p and q be relatively prime and suppose that t is an integer such that 1 < t < 
p + q. Then 

^-) + ^^*g((^))0 + |((^))((5)) 

2pq 2\p q pq J 4 12 \q pq p) \\pJJ \\ Q J J 

□ 

This version of Gessel's reciprocity law bears something surprising: its form is essentially identical 
to Knuth's version [K] of Rademacher's reciprocity law for the Dedekind-Rademacher sums, with 
two differences: Gessel's theorem requires t to be an integer, whereas t E R in Knuth's reciprocity 
law. On the other hand, the conditions in Knuth's theorem are p < q and < t < q, which suffices 
for all practical purposes, however, Gessel's range on t is bigger. In fact, is not hard to remove the 
integrality condition on t, which unifies Knuth's and Gessel's reciprocity theorem: 

Theorem 5.6 Let p and q be relatively prime and suppose that t is a real number such that 1 < 
t < p + q. Then 

f v , v [~t] 2 [~t] (11 1 \ 1 1 

a(q,p,-t)+a(p,q,-t) = L-L + L-l I- + - + — \ + - + — 
v ' v ' 2pq 2 \p q pq J 4 12 

~ ((~P~)) ~ ((~)) ~ {p + q)~4p~~4q" 
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Proof. We will denote the fractional part of x by {x} = x — [x]. Let t = n + r, where n £ Z and 
r £ R, < r < 1 (so r = {t}), then 

Hence 

a(o, 6, t) = a(a, b, [t]) - = a(a, b, [t]) - !((*)) - ^ . 
Now we can use Corollary 15 . 51 for a (a, b, [t]), and the statement follows. □ 
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